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Abstract 

On the basis of the NJL model as an effective theory of QCD and analogies with con- 
densed matter physics, we extract simple physical pictures of the chiral phase transition 
at finite temperature T and/or chemical potential // and predict interesting phenom- 
ena associated with the phase transition which may be seen in relativistic heavy-ion 
collisions to be performed in RHIC or LHC. We indicate that the macroscopic equation 
describing non-equilibrium phenomena of the chiral transition like DCC might be given 
by a diffusion-like equation rather than a wave equation especially when the baryon 
density is present. 



1 Introduction 

The aim of the work to be reported here is to extract simple physical pictures of the chiral 
phase transition at finite temperature T and/or chemical potential fi and to predict inter- 
esting phenomena associated with the phase transition which may be seen in relativistic 
heavy-ion collisions to be performed in RHIC or LHC. The approach adopted will be not 
a fundamental nor very microscopic one like lattice QCD reported by Kanaya|T|; we shall 
make semi-phenomenological arguments based on some effective theory of QCD, the NJL 
model0, and also analogies with possible counterparts of condensed matter physics. A part 
of the present report is based on a review paper ^ based on the NJL model. A remarkable 
point is that the basic picture presented in has been being confirmed in subsequent studies 
with lattice QCD [^] and in particular instanton liquid modelQ. In this report, we shall also 
present further developments and predictions obtained in the approach mentioned above. 

^Talk presented at YITP International Workshop on Physics of Relativistic Heavy Ion Collisions, 9-11 
June 1997, Kyoto, JAPAN. To be pubhshed in Supplement of Prog. Theor. Phys. 
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2 The dependence of thermodynamical quantities on 
the number of active flavors and the vector couphng 



The lattice simulations ^ ^ show that the order and even the existence of the phase transi- 
tion(s) are largely dependent on the number of the active flavors especially when the physical 
current quark masses are used: For rriu ~ rrid ~ lOMeV << lOOMeV ~ rris, the phase tran- 
sition may be weak 1st order or 2nd order or not exist. 

The gross feature of the T dependence and the striking difference between the condensates 
of u (d) quark and the s quark are well described by the NJL model[^ S 0- ^ minimal 
model which describes the low energy phenomena related with the chiral symmetry may be 
given by the generalized Nambu-Jona-Lasinio (NJL) model with the anomaly term]^]; 

C = q[i^ .d~ui)q + J2 ^[{'l>^aqy + (g^A,75g)'] + 9^[detqi{l - ^^)qj + /i.e.], (1) 

a=0 ^ 



where A'^ (a=0~8) are the Gell-Mann matrices with Ao=y |l. SUv{^) symmetry is explicitly 
broken by the current quark masses, m = diag(m„, m^, ms). The determinantal term has 
the SU^i?)) ® SUj^{3) invariance but breaks the Ua{1) symmetry; this term is a reflection 
of the axial anomaly in QCD and induced by instantons[^]. The anomaly term causes a 
mixing in flavors: In the chirally broken phase, it induces effective 4-fermion vertices such as 
{dd){uu){ss) and —{dd){ui'-^^u){si'-y^s), where the former (latter) gives rise to a flavor mixing 
in the scalar (pseudo-scalar) channels. 

The numerical calculation of the thermodynamical potential Q{{{uu)), {{ss))) as a function 
of the quark condensates {{qiqi)) {i = u,s) shows the following||lCI||: In the chiral limit, the 
phase transition is of second order, although the thermodynamical potential is asymmetric 
with respect to the zero condensate due to the cubic term coming from the determinan- 
tal term. With realistic current quark masses {rriu = ma = 5.5 MeV, = 135 MeV), 
Q{{{uu)),0) , i.e., in the u-d sector, has an asymmetric double well, while fi(0, {{ss))), i.e., in 
the strangeness sector, has only single minimum for all temperatures, which implies that the 
phenomenon of the disoriented chiral condensate (DCC)|TT[] could not be expected in the 
strangeness sector in contrast to the non-strange sector. 

It is also noteworthy that the entropy S{T) calculated with the NJL model shows that 
S{T) goes up with with increasing T in accordance with the fact that the chiral restoration is 
a phase transition from an ordered phase to a disordered phase [0; the increase of the entropy 
may imply [0 that the NJL model also qualitatively reproduce the ratios e/T^, (e — 3P)/T^ 



etc obtained in the lattice QCD. 

The temperature dependence of the quark condensates shows that at high temperatures, 
the flavor 5'[/(3)/-symmetry becomes worse 0, which may reflect in the baryon and the 
vector meson spectra, because they are well described by the constituent quark models. 

So much for the case with fi = 0. How about the phase transition at finite density? A 
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numerical calculation [0] shows that at low temperatures lower than about 50 MeV, the phase 
transition is of strong first order in the density direction. Note that our model Lagrangian 
has no vector term like g^i^'^^qY- As a matter of fact, the strength and even the existence 
of the 1st order transition are strongly dependent on the strength of the vector coupling 



W\\ the vector term prevents a high-density state; see also §5. 



3 Mass shift and character change of hadrons 



The NJL model ^ predicts that the mass of the sigma meson decreases as T is raised 
till Tc, a "critical temperature": is found to be constant as long as T < Tc] may be 
defined as the temperature at which starts to go high. This may be in accordance with 
the lattice results on the screening masses [1^. This suggests that at high temperatures the 
decay cr ^ 27r would get suppressed and finally hindered, and then only the electro-magnetic 
process a ^ 27 is allowed. It means that the sigma meson may show up as a sharp resonance 
with the mass ma ~ 2m^. Thus we propose to observe ti~^tt~ , 27r°, 27 and construct the 
invariant mass and examine whether there is a bump in the mass region 300 to 400 MeV. One 
can expect that in the charged system, the process a — 7 — 2 leptons become possible, 
because vr"*" and vr" have different chemical potentials with each other |T8[ . 



The behavior of kaon at finite temperature was examined by the present author using the 
5'?7(3)-NJL model Eq.(l) [0. It was found that as long as the system is in the NG phase, 
the mass of kaon rriKiT) keeps almost a constant, the value at T = 0. 

One can argue that the vector mesons p, uj and <^ decrease their masses as the chiral 
symmetry gets restored; the constituent quark model gives for instance, m,^ ~ 2Ms with Mg 
being the constituent strange quark mass which may be identified with the dynamical one 
generated by the chiral symmetry breaking!^. Other QCD- motivated models and theories 



also predict the decrease of the vector meson masses at finite T |T^, although it is not 
necessarily clear how the decrease of the quark condensates is reflected in the mass shift in 
such theories. 

If the strength of the anomaly term decreases, a character change of r] and 77' mesons is 
expected. A model calculation ||2^ shows that the 77-77' system tends to be well described by 
the non-strange-strange bases (quark bases) rather than the flavor octet-singlet ones at high 
temperatures. This has been confirmed also by instanton liquid model 0. 
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4 Precursory soft modes of chiral transition and chan- 
nel dependence of hadronic correlations at T > Tc 



If a phase transition is of second order or weak first order, there should exist precursory 
soft modes in the symmetric phase prior to the phase transition to a disordered phase; the 
soft modes become tachyonic on the unstable "vacuum" after the critical point. This is well 
known in condensed matter physics and nuclear physics, and there are lots of examples of 



such soft modes; see for an example in nuclear physics. The soft modes are actually 
fluctuations of the order parameter of the phase transition. In the chiral transition, thus 
the soft modes corresponding to the fluctuations {{{ipipY)) and hence {{{ipi'y^Ti/jy)) due to 
the chiral symmetry are expected to exist in the symmetric phase or the high-temperature 



§hase. The NJL model was used to demonstrate explicitly that this is the case [p!4| , |15 
The subsequent lattice simulations on screening masses and studies of the instanton 
molecules support the existence of the soft modes. 

Remarkably, the lattice simulations ^ seemed to show that vector-mesonic and bary- 
onic modes were also obtained in the high-T phase. Actually, the screening masses of the 
vector modes coincide with 2ttT within the error bars, which may simply indicate that the 
interactions between q-q's in this channel are absent or greatly suppressed. The nucleons 
exist as a parity doublet, which was subsequently showed to be compatible with chiral sym- 



metry |23|: A linear sigma model with parity doubling can be constructed which is not in 



contradiction with the low-energy phenomena including the tt-N coupling constant. 

Is there soft modes in the strangeness sector associated with the chiral transition? Soft 
modes are expected to exist when the phase transition is of second order or weak first order, 
so that the restoring force for fluctuations of the order parameter becomes small gradually. 
In the strangeness sector, the thermodynamic potential has single minimum and shows an 
only a shift of the minimum point, which suggests that there is no necessity of a softening 
of the excitation energy for fluctuations of the order in the strangeness sector. A numerical 
calculation shows that this is the case in contrast with the non-strange sector [p^ . 



5 Fate of vector couplings and correlations in the vec- 
tor channel at high temperature 

One may say that there is a reason d'etre of the pionic and sigma-mesonic excitations near 
the transition point even in the high-T phase; they are fluctuations of the order parameter 
of the chiral restoration. Then, is there any fundamental reason why hadronic excitations 
in the vector channel? The approach of the hidden local symmetry p4| seems to claim that 
the existence of the vector mesons are intimately related with the chiral symmetry and its 
spontaneous breaking. We show that the great suppression seen in the screening mass is 

^DeTar also conjectured the existence of such modes in a quite different context [E2[. 
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consistent with the T-dependence of the quark-number susceptibihty Xq{T) obtained by 
the lattice simulations p6[ . 



The quark-number susceptibility Xq is the measure of the response of the quark number 
density to infinitesimal changes in the chemical potentials = u, d) if we confine ourselves 



to the two-flavor case 26, 281; thus 



X,(T, ^) = —p^ = j3j d^iiqiO, x)7og(0, x)g(0, 0)7og(0, 0)ra), (2) 

where Pq = {{3Nb))/V is the quark-number density with Nb being the baryon number 
operator and (3 = 1/T V the volume of the system. We remark that the number susceptibility 
Xq at finite density is directly related to the (iso-thermal) compressibility Kj, as Hj. = Xq/p^- 
Thus one sees that if Xq of a system is large, the system is easy to compress, which may be 
a reflection of a weak repulsion between the constituents of the system. 

The lattice simulations |^ show that as T is raised Xq at /ig = increases very rapidly 



around the critical point of the chiral transition: The ratio of the values in the high-T and 
low-T phase reads Xq(Tc + e, 0)/xq{Tc — e, 0) = 4 ~ 5 with the light two flavors, where e is a 
small number, say 0.03Tc. 

One can easily verify that Xq{T) of the free quark gas increases as M{T) decreases and 
reaches NfT"^ at M{T) = 0: The enhancement is, however, found to be merely about 1.6 
with M{T) as described in the effective theory Eq.(l). Thus there must be an additional 
mechanism to increase Xq to realize the anomalous (relative) enhancement obtained in the 
lattice QCD. One may first note that Xq is the density-density correlation which is nothing 



but the 0-0 component of the vector- vector correlations or fluctuations . Accordingly, Xq 
is further related with the retarded Green's function or the response function in the vector 
channel, the poles of which give the masses of the vector mesons. Thus one recognizes that 
the quark-number susceptibility is intimately related with the properties of the vector mesons 
and fluctuations in the vector channel. 

To demonstrate the relevance of the vector correlations to Xq^ a model calculation with 
a dynamical model was performed [|25| with a simple NJL model having a vector coupling. 
Although one only needs the case with a vanishing chemical potential to compare the results 
with the lattice simulations, it is noteworthy that when fiq there arises a coupling 
between Xq and the scalar-density susceptibility Xs owing to the non-vanishing "vector-scalar 
susceptibility" Xvs^ which are deflned by 

Xs = -^^^|^ = /?/rfx((^(0,x)V^(0,x)V^(0,0)7/;(0,0))), 

dpq 



Xvs = ^^^i^^ = /3/rfx((^^^^ (3) 



respectively. One may also note that Xs is the fluctuation of the order parameter of the 
chiral transition, and is related with the sigma meson propagator. Thus when pq ^ 0, which 
the lattice QCD is difficult to calculate, the properties of the sigma meson reflects in Xq- 
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This is a good example of mode-mode coupling known in the dynamic critical phenomena 
in condensed matter physics. 



We remark that when /i^ = 0, the coupling between Xq ^-^d xs disappear. Putting 
/iq = into the expressions one obtains Xq = X^qHT)/(l + 2gyxfH'^))^ where Xg^H^) is the 
susceptibility for the free-quark gas and Qy the coupling constant in the vector channel as 
mentioned in §2. The denominator of Xq is essentially the inverse of the propagator of the 
vector meson in the ring approximation, at the vanishing four momenta. One sees that Xq 
is suppressed with the vector coupling; is positive. It reasonably implies that the system 
becomes uneasy at /Xg 7^ to compress with the vector coupling; recall that Xq is proportional 
to when /x^ 7^ 0. 



It was shown that the lattice results 26 1 can be accounted for by a decreasing vector 



coupling as well as the decreasing dynamical quark mass due to the chiral transition |^ . Thus 
the lattice results may simply suggest that the vector coupling is suppressed or vanishes in 
the hight-T phase, which implies that the vector correlations and hence the collective vector 
modes would disappear in the chirally restored phase. 

It is interesting that this picture is consistent with the observation that the screening 
masses of the vector modes obtained in the lattice simulations almost coincides with 27rT, 
the lowest screening mass of the q-q system in the chiral limit. 



6 What can we learn about DCC from physics of su- 
perconductors? 



Recently much attention has been paid to time-dependent phenomena of the chiral tran- 
sition, especially in relation to the disoriented chiral condensate (DCC)[^. One of the 
basic questions of this problem is what is a macroscopic phenomenological equation which 
can describe the time-dependent phenomena. Usually, the phenomenological linear sigma 
model is used. This is, however, not necessarily justified. First we note that the problem 
is a dynamical critical phenomena of the chiral transition at finite temperature, and there 
is much resemblance between the chiral transition in QCD and the phase transition to a 
superconductors]^, 0. Abraham and Tsuneto p^] examined what type of equation could 



describe slow and long-wave length phenomena in superconductors; see also ||3^ for sub- 
sequent development. They clarified that the answer depends on the kinematics in which 
one is interested in: For example, in the vicinity of the critical point, the equation is well 
approximated by the non-linear diffusion equation, called time-dependent Ginsburg-Landau 
(TDGL) equation, not a hyperbolic equation as derived from the linear sigma model. The 
origin of the dumping is the Landau damping. It means that it is an open and hence inter- 
esting problem to explore what type of equation is suitable to describe the DCC, and to see 
what kinematical conditions are relevant to relativistic heavy-ion collisions. A preliminary 
calculation shows that when baryon density is present, the relaxation effect due to the Lan- 



dau damping is significant [0] ; we notice that the scalar- vector coupling is also important at 
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finite /i 



7 Summary 

We have seen that the static properties, as represented by the quark condensates {{uu)), {{ss)), 
the thermodynamical potential Q{{{uu)), {{ss))), of the chiral phase transition can be un- 
derstood in terms of a simple model Eq.(l). We have indicated that the flavor dependence 
should reflect in the phenomena of DCC; DCC may be only relevant to the u, d-sector. We 
have also shown that when the baryon density is present, i.e., /i 7^ 0, the strength of the 
vector coupling may change the nature of the chiral transition substantially. 

We have seen that in relation to the mass shift of the sigma meson, the detection of the 
process a —>■ 27t, 27 and 2 is interesting. We have also seen that that rj-r]' mesons may 
change their characters at high temperatures. 

We have also seen that the non-perturbative effects at T > are channel dependent: 
the sigma and pionic channel are special because they are intimately related to the chiral 
symmetry and its spontaneous breaking. On the other hand, the interactions in the vector 
channel are likely to be weak. It was indicated that the parity doubling in the baryon sector 
is consistent with chiral symmetry and the low energy phenomena at zero temperature. 

We have seen that the quark-number susceptibility is intimately related to fluctuations 
in the vector channel in the system: We have suggested that the vector correlations seem 
to be suppressed after the chiral transition. Direct measurements of this behavior of the 
susceptibility in experiment would be extremely interesting. We have also emphasized that 
at finite baryon density, there arises a mode-mode coupling between the scalar and the vector 
channels. 

As for DCC, we have indicated on the basis of an analogy with superconductors that the 
decay process due to the Landau damping is important and it may change the macroscopic 
equation which describes the DCC from the wave-equation like to a diffusion- like equation; 
it also depends on the relevant kinematics. We remark that in non-equilibrium phenomena 
like DCC, effects of the mode-mode coupling will affect the dynamics considerably. 
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